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Abstract 

The Tricritical Ising model perturbed by the subleading energy operator $(3) was 
known to be an Integrable Scattering Theory of massive kinks [14] and in fact preserves 
super symmetry. We consider here the model defined on the half-plane with a boundary 
and computed the associated factorizable boundary S-matrix. The conformal boundary 
conditions of this model were identified and the corresponding S-matrices were found. We 
also show how some of these S-matrices can be perturbed and generate "flows" between 
different boundary conditions. 
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1. Introduction 



Quantum field theory of systems with boundaries have recieved much attention in 
recent years due to substantial advances in several apsects of this subject. One successful 
area of research is in the study of s-wave scattering of electrons from impurities [1] em- 
ploying the technology of boundary conformal field theory [2] . Another fruitful avenue of 
work is the calculation of the boundary S-matrix for the elastic reflection of particles off 
the boundary [3]. In the context of quantum spin chains, this matrix is known as the K 
matrix [4, 5] . Such K matrix facuilties the generalization of the quantum inverse scattering 
method to the case of an open spin chain [6]- [9]. While for the case of integrable perturba- 
tion of conformal field theory, the boundary S-matrix is factorizable if integrable boundary 
conditions are choosen [10]. By applying the technique of thermodynamic Bethe ansatz 
with the boundary S-matirx, many physical results relating to boundary phenomena can 
be obtained [11]- [13]. 

The current work concerns the boundary S-matrix for the tricritical Ising model. 
As is known, perturbing this model in the bulk by the sub-leading energy operator v' 
generates a flow from the c = ^ conformal held theory down to a held theory of massive 
kinks [14]. ^ In fact the resultant theory is supersymmertic and the supercharges can 
be explicitly constructed as integrals of motion [14]. In the present work, we consider 
this model defined on the half-plane with a boundary. The general S-matrix describing 
the factorizable scattering of kinks from the wall is derived by solving the corresponding 
^ The same perturbation with an opposite sign in the coupling constant generates a 
massless flow down to the Ising flxed point [15, 16] 

2 



"boundary Yang-Baxter equation" . The conformal boundary conditions of the model can 
be classified following the proposal of [2]. This is done by using the RSOS picture of the 
Tricritical Ising model. The S-matrix associated with the conformal boundary conditions 
were found and we also studied the perturbation of these boundary conditions by a relevant 
boundary field. 

2. Backround on the Tricritical Ising Model 

The 2-Dimensional Tricritical Ising Model is of interest to physicists because it describes 
tricritical phenomena in a variety of microscopic models [17]. In particular, it is the 
simplest known statistical model to exhibit Supersymmetry [18], and the fact that it can 
be realized experimentally [19] makes it an important model to study. 

Several lattice realization of this tricritical behavior is possible, among them the 
Blume-Capel Quantum Spin Chain [20] and the A4 Restricted Solid-on-Solid (RSOS) 
model [21]. The microscopic model that we will consider here is the Ising model with 
annealed vacancies [22] . The classical Hamiltonian is 



where i,j label the lattice sites, each of the variables ai takes three values —1,0,1 (0 
represents the vacant site) and the first sum in (1) is made over nearest neighbor pairs. 
Concentration of the vacancies is controlled by the chemical potential fi. The phase di- 




(1) 



<i,j> 
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agram of this system in coordinates of temperature T and e~'^ is shown (schematically) 
in Fig.l. There are two phases separated by a line AB of phase transition. For low va- 
cancies and low temperature, we have an ordered phase I with spontaneously broken Z2 
symmetry. On the other side of the transition line we have a disordered phase II with 
unbroken symmetry. The point B is the usual critical point of the Ising model, and the 
solid segment BT corresponds to the second order transition belonging to the Ising model 
universality class. The dotted segment AT is a line of first order phase transition. The 
point T where the two critical segments joined smoothly is the tricritical point of the sys- 
tem. This model belongs to the universality class of the Landau- Ginzburg (/?^-theory at 
this tricritical point [23], so the phase diagram could also be understood in terms of the 
Landau-Ginzburg effective potential of this system [14]. This effective potential (specific 
free energy) is plotted as the function of the order parameter (< a >) of this system in 
Fig.2 for various regions of the phase diagram. The ground state is degenerate in phase I, 
while in phase II it is non-degenerate. There is only one vacuum on the line BT, typical of 
a second order transition. Note that along the first order transition line AT, the vacuum 
is in fact three-fold degenerate. 



At the tricritical point T, the system is described by the conformal field theory (CFT) 
with central charge c = ^ [24]. There are six irreducible representations of the Virasoro 
algebra, and the corresponding conformal dimensions of the primary field (/>(r-,s) are 

organized into a Kac table in Fig. 3. 
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Figure 3 

The fields in this theory can be classified according to their properties under the Z2 
spin- reversal transformation cr — > — cr (this is a symmetry of the microscopic model (1)). 
In particular we have four even fields: the identity / = $(0)5 the leading energy density 
e = ^(^), the subleading energy (vacancy) density e' = ^(^) and the irrelevant field 
e" = ^(|)- Here we use the short-hand notation for the scalar field ^(a,a) with 

A = A. There are two odd fields under spin- reversal: the leading spin field o = ^(^): and 
the subleading spin field a' = ^(X)- 

In addition to the scalar fields, this conformal field theory also possesses fermion 
fields G{z) = ^(3,0) ^(^) = ^(o,|)- The currents G{z) {G{z)) together with the 
stress tensor T{z) {T{z)) generate the Neveu-Schwartz-Ramond algebra, hence this model 
exhibits superconformal symmetry [18]. We can alternatively classify the fields in this CFT 
according to this extended symmetry. For example, in the Neveu-Schwartz sector of this 
CFT we have the superfield 

where = _^-), and = |) are fermion and antifermion fields respectively. All 
the component fields of ^(^) are mutually local and, together with the identity operator 
/ and its descendants, they constitute the Neveu-Schwartz sector of this CFT. In this 
sector, Kramers- Wannier duality transformation acts as a second Z2 symmetry under 
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which *(^) — > — *(^); — > —9. More specifically, under duality transformation, the 
signs of e, and G are changed, while /, e', and G are unaltered. 

The spin fields a and cr' provide the Ramond representation of the superconformal 
symmetry [18]. The duality transformation converts them into the disorder fields ji and 
jj! respectively [25], where /x = Gqg^ ji' — Goa' . Here Go is the zero-mode of G in the 
Ramond sector. Hence the Operator Product Algebras generated by (/, cr, cr', e, e') and 
{I, jj,, jj,' , —e, e') are isomorphic under duality. In other words, the Trier itical Ising Model 
is a self-dual system. 



3. Supersymmetric perturbation of the Tricritical Ising Model 



In the c= jq conformal field theory, the energy densities and spin fields are relevant 
operators. Perturbation of the bulk theory by these operators were studied numerically 
in [26]. This work confirms the conjecture that perturbations by the operators e, e' and 
cr' are individually integrable [27]. In particular, the transition line AB in Fig.l can be 
generated by perturbing the tricritical fixed point Hamiltonian by the vacancy density e' 
[15]. Thus the scaling region around the point T can be described by the field theory with 
the "action" 

H^ = H^^^^+X J ^^i)d'x, (2) 

where H^j_^ is the "action" of the bulk CFT. The renormalization group (RG) flow from 
T down to the Ising fixed point B is given by (2) with A > and can be described by a 
massless scattering field theory [16]. For A < 0, (2) is an integrable massive field theory 

6 



[14] and is associated with the first order transition fine AT. According to [27], this model 
possesses an infinite number of bosonic Integrals of Motion, and we shall concentrate on 
the case with A < 0. Further more, because $(3) is the upper component of the superfield 
^(_^), so the perturbation (5.2) does not destroy the global supersymmetry of H(^7_y In 
fact the supercharges 

Q = J [Gdz + AXi)dz\; Q = J [Gdz + AXipdz] 

are fractional spin integrals of motion in the perturbed theory [14] . They were used in [14] 
to construct the factorizable scattering S-matrix given below. 

Along the line AT, the vacuum is three-fold degenerate, thus the particle spectrum of 
(2) with A < must contain "kinks" separating the domains of different vacua. As shown 
in Fig. 4a, if we label the vacua as —1, 0, 1, then there are four sorts of kinks of mass m: 
Kq^+i, Kq^_i, Ki^o, K_ifi. An asymptotic N-kink scattering state, say an in-state 

|-^CTo,ai(^l)-f^ai,a2(^2)-f^a2,a3(^3)---f^cTjv-i,CTjv(^Af) >in (3) 

where 6i is the rapidity of the z-th kink, can be thought of as the sequence of N + 1 
domains of vacua uo, (Ji, (Jat placed along the x-axis, with the kink i^o-i.o-i+i separating 
two neighboring domains and CTj+i. The neighboring vacua must satisfy the "admis- 
sibility condition" |crj — crj+i| = 1. Thus in (3) the vacua +1 and —1 must be separated 
by the vacuum 0. The out-state of N-kink scattering is related to the in-stated (3) by the 
factorizable S-matrix built up from the two kink S-matrix. This S-matrix consists of four 
different amplitudes (shown in Fig. 4b) Ao{9) , Ai{9) , Bq{9) , Bi{9) , defined as 

\KoAdl)KsA(^2) >.n= Ao(^l2)|i^O,.(^2)i^.,o(^l) >out +Ai (^12) |i^0,-. (^2)i^-«,0 (^l) >out; 

(4a) 



|i^.,0(^l)i^0,.(^2) >in= Bo{9l2)\Ks,o{e2)Ko,s{ei) >ouU W 
\Ks,o{ei)Ko,-s{e2) >in= Biiei2)\K,,o{e2)Ko,_M >ouU (4c) 

where s = +1 or —1, and 6\2 = Oi — ^2- These amphtudes are given in [14] as 

9 9 
Aq = cosh{-)A{9); Ai = -isinh{-)A{9); 

Bo = cosh{^--'-^)Biey, B, = cosh{^- + '-^)Bie). (5) 
The "minimal" expressions for A{9) and B{9) are 

A{e) = exp{^9)S{ey, B{9) = V2exp{--f9)S{e), (6) 

where exp{27rij) = 2 and 

Note that these amphtudes do not exhibit any poles in the "physical strip", hence there 
are no bounded state in the bulk scattering theory. 

4. The General Boundary S-Matrix 

In this section we consider the perturbed Tricritical Ising model in the semi-infinite 
plane {{x,y) : x < 0; — oo < y < oo}. The general symbolic action with general boundary 
condition has the form 

/oo pO pco 

dy ^z{x,y)dx+ ^B{y)dy, (8). 
-oo J —oo J —oo 



where H7__^qqq is the CFT action on the half-space with Conformal Boundary Conditions 
(CBC) on the boundary x — Here we assume a choice for the boundary operator $b 
can be made such that the integrabihty of the bulk theory is preserved. Therefore the field 
theory (8) can be described by a factorizable scattering theory of kinks with the boundary. 
Since the bulk theory contains three degenerate vacua, we can expect the boundary to 
exist in one or more of these vacuum states. In the general scenario, the boundary states 
need not be degenerate, and some of these states can appear as boundary bound states 
in the scattering process [10]. If some of the boundary vacua are degenerate, then the 
scattering of kinks with the boundary can change its state. 

Let us denote the "boundary creating operator" by Ba where a G {—1, 0, 1} labels the 
boundary vacuum state. The Fock space of the boundary scattering theory consists of the 
scattering in-states 

\K^,,a,{ei)K^,,^,{92)...K^^,aBa >^n (9) 

with 6*1 > ^2 > ••• > ^AT > 0. The out-state of the scattering is again related to the in-state 
(9) by a S-matrix built up from the two-kink S-matrix (4) and the one-kink boundary scat- 
tering S-matrix. The boundary S-matrix comprised of six elementary amplitudes defined 
as 

\K+,,oid)Bo >in= R+mK+i,oi-^)Bo >ouu (10a) 
|K_i,o(^)i?o >m= R-mK-iA-^)Bo >ouU (106) 

|Ko,+i(^)s+i >in= p+{e)\Ko,+i{-e)B+i >out +v+{e)\Ko,-i{-e)B_i >out; (loc) 

\Ko,-iid)B_i >in= P-{d)\Ko,-ii-e)B_i >out +V-{e)\Ko,+ii-d)B+i >out, (lOrf) 
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and the associated space-time diagrams are illustrated in Fig. 5. Note that in the general 
case we do not assume a priori that the boundary will respect the spin-reversal symmetry 
(cr — > —a) of the bulk theory. 

Factorizability of the scattering theory (8) means that the boundary amplitudes (10) 
have to satisfy the Boundary Yang-Baxter equation [3, 10] (shown in Fig. 6) and results in 
the following set of functional equations: 



R+AoR+Ai + R+AiR-Ao = R-AiR+Aq + R-AqR-Ai; (11a) 

P+BoV+B^ + V+B^P_Bi = V+BoP+Bo + P_B^V+Bo; (116) 

P+BiV+Bo + V+BoP_Bo = V+B^P+B^ + P-BoV+Bi; (11c) 

P+BiP+Bi + V+BqV-Bi = P+BiP+Bi + V-BqV+Bi; {lid) 

P+BoP+Bo + V+B^V_Bo = P+BoP+Bo + V_B^V+Bo, (lie) 



(the other equations are obtained by interchanging -|- — in (11)) where each term in the 
equations has the arguments Ri{9i)Sj{92 + 9i)Rk{02)Si{92 — 9i). Notice that the equation 
for and R- decouples from the equations for the other amplitudes. This is to be 
expected since the scattering associated with R+ and R- does not involve changes in the 
boundary state. However if the amplitudes (11) exhibit "boundary bound states", then 
the "boundary bound-state bootstrap equation" will mix the various amplitudes. This 
will be discussed in Section 6. Solution to these equations can be obtained by elementary 
methods and has the form 

R+{9) = {l + Ash^)M{9y, R_{9) = {1- Ash^)M{9); (12a) 
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P+{9) = {X + Ysh9)N{9y, P_{9) = {X -Ysh9)N{9); (126) 

V+{9) = k+sh^N{9); V_{9) = k_sh^N(9), (12c) 

where A, X, Y, k+ and A;_ are constants depending on the boundary condition. It is worth 
noting that R+ and i?_ contains one free parameter A, while the other amplitudes involve 
three independent parameters (we can always absorb a constant into N{9) when needed). 
The normalization functions N{9) and M{9) can be determined by the boundary unitarity 
and cross-unitarity constraint for the boundary S-matrix. For this model, the unitarity 
conditions assume the form 

R+(9)R+{-9) = 1; (13a) 

P+{9)V+{-9) + V+{9)P-{-9) = 0; (136) 

P+{9)P+{-9) + V+{9)V-{-9) = 1; (13c) 

and other equations with the substitution + -s-^ — in (13). 

To obtain the cross-unitarity condition, we have to analytically continue the boundary 
S-matrix (10) into the cross channel [10], resulting in the following equations 

-0) = Aoi29)R+{^ +9)+ ^i(2^)i?_(^ + 9); (14a) 

^'+(y-^) = 5o(2^)P+(y+^); (146) 

V+{--9)^B^i29)V+{- + 9), (14c) 

and similar equations under the interchange + <-> — in (14). The unitarity conditions (13) 
and (14) reduce to 

mO)M{-9)= (15a) 
11 



(156) 



N{e)N{-e) 



1 



(15c) 



-Y'^sh^e-k+k-sh'^^' 



N{--9)^Bo{29)N{- + 9), 



(15d) 



for the normalization functions N{9) and M{9). The solution to these equations re- 
quires analytic information about the scattering amplitudes (12) which are governed by 
the boundary condition. More specifically we need to determine the existence, if any, of 
"physical poles" (ie poles in the physical strip) in these amplitudes. To this end we will 
describe in the next sections the various conformal boundary conditions in this model and 
their associated scattering S-matrices. 

5. Conformal Boundary Conditions of c = ^ CFT 

A vareity of boundary conditions are possible for a given bulk theory, and they can 
be catagolrized into universality classes much like the case of the bulk theory. Each uni- 
versality class of boundary conditions are characterized by a boundary fixed point [28] 
which is invariant under conformal transformations. Various conformal fixed points can be 
connected by renormalization group trajectories. The conformal boundary conditions for 
c= jQ may be classified following the work in [2] . Conformal invariance of the boundary 
is equivalent to the requirement that the two components T and T of the stress tensor be 
equal on the boundary. If we denote a conformal boundary state by |a), then it must be 
annihilated by (L„ — L„) for each n. Solutions to this constraint were given in [29], and 
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they are linear combinations of the states 

|i>=5^|i,iV>®|j,iV>, (16) 



AT 



where j labels a highest weight representation of the algebra of the and \j,N > is an 
orthonormal basis in this representation space. Similarly, \j, N > forms an orthonormal 
basis in the j representation of the algebra of L„. ^ The boundary states corresponding 
to "physical" boundary conditions can be constructed following the procedure in [2] and 
have the form 

|0) = C[\o) + v\^) + v\l) + If) + ^1^) + ^li)]; 

li) = C[rj'\0)-V-'\^o)-V-'\l)+v'\l) - ^V'l^e) + ^^"^1^)]; 
I) = C[v'\0) - + r/^|i) + ^r/^l^) - ^ry-^l^)]; 

I) = Cm+r^\^)+f^\i) + If) - ^1^) - ^li)]; 
4) = V2Cm - v\^o)+v\i) - If)] = ^,m\0) + If))]; 

II) = V2C[v'\0) + v-'\^o) - - = j^mi) + If))]' (17) 
where C = ^ and r] = ^ ^sin^ • Here V denotes the Kramers- Wannier duality 
transform. Thus we can regard the state |^) as the dual boundary condition of |0) and 
||). Similarly |^) is dual to ||) and \jq)- 

To interpret the boundary conditions associated with the states in (17), it will be 
illustrative to refer to the A4 RSOS lattice realization of this model [21]. Consider a 
diagonal square lattice where the degree of freedom at each site can take on four values 



^ Alternatively, we could try to classified the boundary states using the super currents G 
and G of the bulk theory [2] . The resultant states will carry highest weight representations 
of the superconformal symmetry. We will not pursue this line of analysis here. 
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li = 1, 2, 3 and 4. The value of li is further constrained by the requirement |Zj — lj\ = l for 
nearest neighbor sites i and j. In addition there is a ferromagnetic next-nearest- neighbor 
interaction and a single-site interaction. Thus the lattice divides into two sublattices I 
and II, where the values of li are odd in one sublattice, and even in the other sublattice. 
It is clear that there are three degenerate ground states as shown in Fig. 9. To make 
contact with the order parameter < cr > of the Ising model with vacancies (1), we shall 
label the ground states in Fig. 9 by —1, and The Hamiltonian of the RSOS model is 
symmetric under the global Z2 transform li — > 5 — li, which corresponds to the spin- reversal 
symmetry of the model (1). The three ground states will become identical in the bulk at 
the tricritical point. However the nature of the order parameter at the boundary will 
depend on the specific boundary condition. As shown in [2], the state |A(i ,5)) corresponds 
to the boundary condition where all boundary degrees of freedom are fixed to the value s, 
this is illustrated in Fig. 10. For the state |0), the boundary degrees of freedom are fixed 
to 1, and hence the neighboring state must be in the state 2. Thus in the continuum 
limit, the order parameter at the boundary will be in the —1 vacuum, and we shall denote 
this boundary condition by (— ). The boundary state |^) corresponds to fixing the the 
boundary degree of freedom to 2, and the neighboring sites can be in either states 1 or 3. 
In this case the order parameter at the boundary may be in the —1 or vacuum. This 
boundary condition where the —1 and vacua are degenerate will be labeled (—0). We 
can similarly associated ||) with the boundary condition (0+) where the and -|-1 vacua 
are degenerate at the boundary. In the same way, the boundary condition {+) where the 
order parameter is fixed to the +1 vacuum will correspond to the boundary state ||). One 
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can regard (— ) and (+) as "fixed" boundary conditions for the Ising model with vacancies 
(!)• 

In the work of [30], they identified the boundary states jA^^ i) > as the boundary 
condition where the boundary degrees of freedom are fixed to the value r, while the neigh- 
boring sites must be in the state r + 1 (for 1 < r < 3). For the boundary state |^), this 
means that the boundary will be in state 2 and the neighboring sites will be in state 3. This 
of course fixes the order parameter to be in the vacuum in the continuum limit, and we 
denote this the (0) boundary condition. A few comments about this boundary condition 
is in order. Using the states in (17), the partition function of the model on the cylinder 
with boundaries at both ends can be computed. Let us define the modular parameter as 
q = e^'^*'^, where r = on a cylinder of length L and circumference R. The partition 
functions involving the boundary condition (0) include 

^(o)(o)(?) = Xo(?) +X|(9), (18a) 

and 

Z(_)(o)(9)=xx(g). (186) 

The result (18a) was conjectured earlier in [31] for the Tricritical Ising partition function 
with microscopic "free" boundary spins, and later confirmed by direct calculation [30] and 
numerically in [32, 33]. Furthermore (18b) was also obtained in [33] by considering the 
partition function with "mixed" boundary condition, ie microscopic "fixed" and "free" 
boundary conditions. Thus we concluded that the boundary condition (0) is associated 
with the lattice Tricritical Ising model where the microscopic boundary spins are not 
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constrained. Such lattice boundary conditions are often referred to in the literature as the 
"free" boundary condition, even in the conformal field theory. 

Lastly we consider the conformal boundary state | ^ ) which has not previously been 
identified. Note that it is dual to the boundary conditions (—0) and (0+), and we conjecture 
that it corresponds to the case where the boundary can exist in all three vacua. In other 
words, like in the bulk theory, the three vacua are also degenerate on the boundary. We 
shall call such boundary condition "degenerate" and label it as (d). Since a great deal of 
"fine-tuning" of the boundary parameters are required to achieve three-fold degeneracy of 
the vacua, we expect the boundary condition (d) to be unstable under perturbation by 
relevant boundary operators. This boundary condition carries the ^ representation of the 
Virasoro algebra, so the boundary operators that can appear along such a boundary is 
given by the Operator Product Expansion of the primary field $(^) with itself [18]. Prom 
the operator product 

we see that two relevant boundary fields •^(-s-) and i^^^) can appear for this boundary 
condition. Perturbing by these relevant operators could generate renormalization group 
fiow from the boundary condition (d) down to a more "stable" boundary fixed point. To 
lend support for this speculation, we computed the boundary ground state degeneracy (the 
so-called "^r-f actor" [34]) from the boundary states (17) and they are in order of decreasing 
magnitude 

9{d) = V2r]'^C; g^_o) = g^o+) = rfC, 

9{0) = V2C; and 5r(_) = g(+) = C. 
16 



The value of g should decrease along renormalization group trajectories, and in this re- 
spect, (d) would be the most "unstable" conformal boundary condition. Similar comments 
can be made about the two-fold degenerate boundary conditions (—0) and (0+) with the 
next-largest value of g. One can repeat the analysis for the associated boundary operator 
for these boundary conditions, and we found that they both admit the relevant boundary 
operator V'd)- Thus we expect that the perturbation of (—0) and (0+) by this operator 
would generate RG flow down to other conformal boundary conditions. This is studied in 
the next section. 

6. Boundary S-matrix for the Conformal Boundary Conditions 

Boundary conditions can in general be changed by altering the coupling constants 
of model on the boundary. A simple example would be the introduction of a "magnetic 
held" which favors a particular direction of lattice spins. In the continuum limit, this 
magnetic field would shift the specific free energy of one of the vacua and destroy the 
ground state degeneracy. Let us consider the massive Tricritical Ising model with the 
conformal boundary condition (—0) and perturbed by the relevant boundary operator 
V'd-). The symbolic action for the corresponding field theory is (8) with 

Hr_+cBC = ^^"^+(-0); = /i'i/'d), (19) 

and h is some boundary coupling constant with dimension [length]~ i . Since the boundary 
operator has the same dimension as the bulk perturbing operator ^(|), this boundary 
perturbation is integrable [10]. The boundary field ip^s^ couples to the boundary spins for 
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the boundary condition (—0) and would alter the specific free energy of the —1 vacuum. 
Thus for h 0, the ground state will no longer be degenerate. 

To be more specific, let us treat the case of a RG flow from (—0) down to the "fixed" 
boundary condition (— ). For h > we expect the vacuum —1 to be the ground state of 
the boundary, while the vacuum is an excited state. The associated boundary scattering 
S- matrix will involve the amplitude P-{$), and since the boundary cannot exist in the 
vacuum +1, we have A;_ = A;+ = in this case. The Yang-Baxter equation (11) becomes 
an identity for this boundary condition, and the amplitude P-{9) will be determined by 
the unitarity conditions (13) and (14). The solution can be written in the form 

P_(^) = Pf^^(^)P^,,(^), (20) 

where Pminifi) is the minimal solution of the equations 

Pmin{0)Pmin{—9) = 1; (21a) 



77r 77r 
Pmini^ - ^) = So(2^)P™n(y + 0), (216) 

with no poles in the physical strip, and has the expression 

'"'"^'^ - L\ r(. + ^)^r(. + i - J^)V(U - 1 - ^) 



The other factor 



- sin^ + ishO' ^^^^ 



is a CDD factor [35] which exhibits a pole at 6* = i^. When this pole lies inside the 
"physical strip" , < ^ < , it could be interpreted as a "boundary bound state" as shown 
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in Fig.ll. If we denote Co to be the specific energy (per unit boundary length) of the 
boundary state \Ba >, then 

Co — e_i = mcos^, (24) 

defines the binding energy of the bound state. Clearly the parameter ^ is related to the 
boundary coupling constant h, and we can analyzes the effect of the RG flow by varying 
the value of ^. Furthermore we can consider scattering of this bound state with other kinks 
by the "boundary bound state bootstrap equation" [10] ^ 

R+{9) ^Bi{9- iOP-{0)Bi{9 + iO; (25a) 

and 

R.{9)^Bo{9-iOP-{9)Bo{9 + iO, (256) 

as illustrated in Fig. 12. This gives two more scattering amplitudes which can occur when 
the boundary is in the vacuum state 0. Equation (25) can be simplified to a form resembling 
the general solution (12): 

R+{9) = ^{cos^+ish^)B{9-iOB{9 + iOP^''''{9)Pmin{9); (26a) 

and 

R_{9) = l{cos^-ish^)B{9-iOB{9 + iOP^''''{9)Pmin{e). (266) 

For the boundary condition (—0) the vacua —1 and are degenerate on the boundary and 
we expect the boundary to admit a bound state at ^ = ^. Thus the amplitudes (20) and 
^ Similar "fusion" procedure was used in [36] to derive the boundary S-matrix for the 
Sine-Gordon model at the supersymmetric points 
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(26) with C = f constitute the boundary S- matrix for this unperturbed boundary condition 
(— 0). When we "turn on" the boundary perturbation {h > 0), (cq — e_i) > and the 
amphtude P-{9) wiU exhibit a pole in the physical strip at ^ = for < ^ < One can 
think of the bound state as formally "gluing" a kink Kq_i of rapidity (interpreted as a 
fluctuating domain wall) to the boundary \B-i >. As we increase the perturbation, ^ ^ 0, 
the bound state becomes weakly bound and the domain wall develops large fluctuations 
which propagate well into the bulk. At ^ = 0, the domain wall has zero rapidity and is no 
longer bounded to the boundary. 

For — f < ^ < 0, the pole leaves the physical strip. In this regime (cq — e_i) > m, so 
there is no "boundary bound state" and —1 remains as the only stable ground state of the 
boundary. At ^ = — |, the poles and zeros of P^^^{9) become two- fold degenerate. As we 
increase h further, ^ develops an imaginary component: ^ = + I'd and the amplitude 
P-{0) exhibits two poles at ^ = — ± These poles depart to inflnity when h gets 
increasingly large. In this limit, the energy difference between the two vacua (eo — e_i) 
becomes inflnitely large and we expect the resultant boundary condition to be the "flxed" 
case (— ). Thus the amplitude P_ will become the boundary S-matrix for the (— ) boundary 
condition in this limit, and we found 

R^_){e) = Pmin{0). (27) 

It is not surprising that this "flow" from a two-fold degenerate boundary condition down 
to a unique ground state resembles closely the "flow" from free to flxed boundary condition 
in the Ising model [10]. 

It is interesting to consider the perturbed action (8) for /j, < 0. In this case we expect 
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the perturbation to "raise" the specific boundary energy of the —1 vacuum, and the vacuum 
would now become the unique ground state. This situation corresponds to the analytic 
continuation of ^ above the degenerate point ^. In the regime ^ < ^ < tt, the difference 
in boundary energies become 

(e_i — Co) = mcos( > 0, (28) 

where C — — Hence we can regard the —1 vacuum as a boundary excited state and 
it could appear as a boundary bound state in the scattering S-matrix. Indeed one can 
check that the amplitude R-{0) now possesses a bound state pole in the physical strip 
at ^ = iC,. The other amplitude R+{0) does not have this pole, consistent with the fact 
that the boundary cannot be in the vacuum state +1. The analysis in this regime mirrors 
that of the domain < ^ < ^, except that the role of —1 and are interchanged. At 
^ = TT, the bound state becomes unbounded and emerges as the only stable boundary 
ground state. As we decrease the value of h further such that ^ > tt, R+{9) and R-{9) 
are the only physical boundary scattering amplitudes and there is no more bound state 
pole. When the RG flow reaches the point for which ^ = 27r, the poles and zeros in the 
factor (cos| ± ish^)B{e + i)B{e - O^f ^^(^) (common to both and R-) become 
two-fold degenerate. Further perturbation would introduce an imaginary component to 
^: ^ = 211 -\- iQ. The amplitudes R+{0) and R-{9) would correspondingly exhibit two 
poles at ^ = 27rz ± g. These poles will separate and depart to infinity when h approach 
negative infinity. Under this RG flow, we expect the boundary specific energy e_i to 
become increasingly large, and thus the boundary will be "fixed" to the stable ground 
state 0. This is the (0) boundary condition and in this limit R+ and i?_ will both be 
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reduced to the boundary S-matrix for (0): 

R^o){0) = e-'"'^Pmin{0). (29) 

The RG flows from (—0) boundary condition down to (— ) and (0) are shown schematicaUy 
in Fig. 13. Of course we could also consider the perturbation of the (0+) boundary condition 
by '03 and we expect it to generate the same flow pattern and conformal boundary S- 

5 

matrices. It is not easy to relate the perturbation coupling h to the S-matrix parameter ^. 
However, based upon the pattern of flow of the scattering amplitudes, we conjecture the 
relation in this case to be 

2, ,2^ TV. 27r, , , 

n. ~ ms [cos(— + — ) — cos— J. (30) 



As shown in Fig. 13, by applying the duality transformation to the RG flows from 
(— 0), we obtained the RG flows from the (d) boundary condition down to either (0) or 
a superposition of (— ) and (+). Note that the held does not change sign under 
duality. The symbolic action of this dual theory is (8), where H7__^QQfj = H^_^(^^^ is 
the CFT action with conformal boundary condition (d). For this new RG flow, the point 
h — corresponds to the boundary condition {d) where all three vacua are degenerate. 
This boundary condition respects the full symmetry of the bulk theory, so the associated 
boundary S-matrix should be invariant under spin-reversal. Thus we have A = Y = and 
k+ = in the general solution (12), and denotes 

R+ = R- = R; P+ = P_ = P and V+ = V- = V. 

When this conformal boundary condition is perturbed with h > 0, the boundary energies 
e+ and e_ are both higher than cq, and hence the vacua —1 and -|-1 are the excited states 
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of the boundary. Here we shall assume the perturbation does not destroy the Zi spin- 
reversal symmetry of (d). The amplitude R{&) can be obtained by solving the unitarity 
constraints (13) and (14). It must exhibit a bound state pole in the physical strip, say at 
iv^ corresponding to these excited boundary states, and has the form 

R{Q) = e-2^^Pf ^^(^)P^,„(^); (31) 

where P^^^{0) is the CDD factor (23). The other amplitudes are again related to R{9) 
by the "boundary bound state bootstrap equation" 

P{e) = [Ao{d - v)Aoie + v) + Ai{e- v)Ai{d + v)]R{e), (32a) 

and 

v{e) = [Ao{d - v)Ai{e + v) + Ai{e- v)Ao{e + v)]R{e). {32b) 

The expressions for P and V can be simplified to 

P{9) = ch^e-^''^A{e - v)A{e + v)P^''''{e)Pmin{0); (33a) 

V{e) = -ish^-e-^-^'Aie - v)A{e + v)P^''''{e)Pmin{0), (336) 

which are the general solution (12) with X — ch^ and — k- — —i. 

For the boundary condition (d), the three vacua are degenerate and all scattering 
amplitudes should possess a pole at ^ = The associated boundary S- matrix is given 
by (31) and (33) with v = ^. Perturbing this boundary condition in (8) with h > would 
lower Co, while the other vacua —1 and +1 will appear as bound states in the amplitude 
-R at ^ = iv. In the domain — ^ < t; < 0, there is no bound state, and R remains as 
the only physical amplitude. Under further RG flow, v develops an imaginary component 
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{v = + ivo) and travels to infinity. One can check that in this hmit, R reduces to the 
amphtude i?(o) ior the boundary condition (0). 

Perturbing (d) with h < raises the energy eo above e+ and e_. For ^ < v < n the 
vacuum enters as the bound state in P and V. In the regime tt < v < 27r, the bound state 
leaves the physical strip, with +1 and —1 remaining as the only stable boundary ground 
states. At ^; = 271, the poles and zeros of the t;-dependent factors of both P and V become 
two-fold degenerate. Further RG flow would introduce an ever increasing imaginary part 
to V, ie V = 2iT + iv' , and drives it to infinity. In this limit, the amplitude V approaches 
zero, while P flows to the "fixed" boundary condition amplitude (27). This means that 
the scattering theory decouples into two equivalent sectors, and the boundary becomes a 
superposition of the pure boundary states |-B+i > and >. Since this RG fiow pattern 
is the same as the fiow from the (—0) boundary condition, we expect the parameter v to 
be related to the coupling h in the same manner as (30). 

7. Discussion 

In this work we computed the boundary S-matrix associated with the conformal 
boundary conditions of the Tricritical Ising model. The interpolating S-matrix which 
describe the renormalization group fiow between various boundary conditions was also 
found. However we have not yet found an integrable boundary perturbation which give 
rise to the general S-matrix (12). It would be interesting to determine this more general 
perturbation. Since the general solution (12) contains more than one boundary param- 
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eter, it is conceivable that the associated perturbation involves more than one boundary- 
operator, and some of them would break the spin-reversal symmetry. Another interesting 
follow-up to this work would be to consider the massless limits of the S-matrix. Since the 
interpolating S-matrix contains a boundary parameter, its massless limit might be non- 
trivial. This information would be relevant to the interpretation of the c — jq conformal 
field theory as a massless scattering theory [16, 37]. 

Finally, we have not discussed in this work the possibility of a supersymmetric bound- 
ary condition [36] . Since the degenerate boundary condition (d) respects the full symmetry 
of the bulk theory, it is plausible that it might also preserve supersymmetry. One would 
like to check if any other conformal boundary states (or a linear combination of them) 
could also be supersymmetric. Furthermore, if some of the conformal boundary states do 
possess supersymmetry, it would be interesting to see if it could survive under perturba- 
tion, for example, by the boundary field V'd)- the level of the boundary scattering 
theory, one could check supersymmetry by determining the appropriate linear combina- 
tion of supercharges Q and Q which survive as an integral of motion. This combination 
should commute with the supersymmetric boundary S-matrix, and could possibly involve 
the topological charge of the massive theory [36] . 
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Figure Captions 

Fig.5.1. Phase diagram of the Tricritical Ising model. 

Fig. 5. 2. Landau-Ginzburg effective potential for the various regions of the phase 
diagram. 

Fig. 5. 3. Kac table of conformal dimensions for c = 

Fig. 5. 4. Space-time diagrams of the four variety of the kinks (a) and their bulk 
scattering amplitudes (b). 

Fig. 5. 5. Space-time diagrams for the General Boundary Scattering Amplitudes. 

Fig. 5. 6. Figure of the "Boundary Yang-Baxter" equation satisfied by the boundary 
S-matrix. Here a,6, c e —1}. 

Fig. 5. 7. Space-time diagram for the Bundary Unitarity equation with a,d & {+1; 

Fig. 5. 8. Figues of the Cross-Unitarity condition with a, 6 G —1}. 

Fig.5.9. The three ground states of the RSOS model and the associated vacuum states 
in the continuum limit. 

Fig.5.10. Various conformal conditions of the Tricritical Ising model in terms of the 
RSOS picture. 

Fig. 5.11. Bound state pole exhibit by the amplitude P_. 

Fig.5.12. Figure of the "Boundary Bound State Bootstrap" equation with 

a e {+1,-1}. 

Fig. 5. 13. Schematic diagram for the renormalization group flow from the (—0) and 
(d) boundary conditions. Here V denotes duality transformation. 

Fig.5.14. Figure of the "Boundary Bound State Bootstrap" equation for the pertur- 
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bation of the degenerate boundary condition. Here a,b E {+1, — !}• 
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Fig. 4 b 
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